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We perform a detailed analysis of the latest CMB measurements (including BOOMERaNG, DASI, 
Maxima and CBI), both alone and jointly with other cosmological data sets involving, e.g., galaxy 
clustering and the Lyman Alpha Forest. We first address the question of whether the CMB data 
are internally consistent once calibration and beam uncertainties are taken into account, performing 
a series of statistical tests. With a few minor caveats, our answer is yes, and we compress all data 
into a single set of 24 bandpowers with associated covariance matrix and window functions. 
We then compute joint constraints on the 11 parameters of the "standard" adiabatic inflationary 
cosmological model. Out best fit model passes a series of physical consistency checks and agrees with 
essentially all currently available cosmological data. In addition to sharp constraints on the cosmic 
matter budget in good agreement with those of the BOOMERaNG, DASI and Maxima teams, we 
obtain a heaviest neutrino mass range 0.04 — 4.2 eV and the sharpest constraints to date on gravity 
waves which (together with preference for a slight red-tilt) favors "small-field" inflation models. 



I. INTRODUCTION 

The recent discovery |l]-||] of multiple peaks in the 
cosmic microwave background (CMB) power spectrum 
marks a major milestone in cosmology. Confirming 1970 
predictions of Peebles & Yu and Sunyaev & Zeldovich 
H , it greatly bolsters the credibility of the emerging stan- 
dard model of cosmology, and allows many of its free 
parameters to be measured with a precision that cosmol- 
ogists have yet to get accustomed to. 

This new precision also offers new opportunities for 
consistency tests, both for systematic errors that might 
affect individual measurements and for incorrect assump- 
tions about the underlying physical processes. The goal 
of the present paper is to carry out these two types of 
tests. 

We begin in Section [il] by testing measurements of the 
CMB power spectrum for consistency with a series of 
quantitative statistical tests, including the effects of cal- 
ibration and beam uncertainties. Since the customary 
plot of available measurements has now evolved into a 
bewildering zoo of over 100 band power measurements 
which is increasingly difficult to interpret visually be- 
cause of calibration and beam effects, we perform an 
essentially lossless data compression of all data into a 
single set of 24 bandpowers with associated covariance 
matrix and window functions, useful as a starting point 
for parameter fitting. 

In Section [II, we compute quantitative constraints on 
the 11 parameters of the "standard" adiabatic inflation- 
ary cosmological model in a variety of different ways, 
using data from, e.g., the CMB, galaxy clustering, the 
Lyman Alpha Forest, Big Bang nucleosynthesis (BBN) 
and Hubble constant measurements in various combina- 
tions. This enables us to identify a number of parameter 
constraints that are robust and consistent with all data, 
as well as areas where there is slight tension between data 
sets pulling in different directions. Although numerous 



such studies have been performed in the recent literature, 
e.g., [@-^|, the dramatically improved precision allowed 
by new BOOMERaNG §, DASI @, Maxima @ and 
CBI p(| data makes it worthwhile and timely to revisit 
this issueQ. The present work extends the solid recent 
analyses of the experimental teams ]l|, p9| , p0| mainly in 
the following ways: 

1. Inclusion of more parameters allows us to place 
constraints on neutrinos and gravity waves and to 
quantify the additional degeneracies that they in- 
troduce. 

2. Joint analysis of all CMB data allows us to place 
stronger constraints and perform consistency tests. 

3. Inclusion of explicit power spectrum modeling for 
the galaxy clustering and for the Lyman Alpha For- 
est allows stronger constraints and important new 
cross-checks. 



II. IS THE CMB STORY CONSISTENT? 
COMPARING AND COMBINING 
MEASUREMENTS OF THE ANGULAR POWER 
SPECTRUM 

In this section we test the CMB data for internal con- 
sistency and combine them into a single set of band pow- 
ers, calibrating the experiments against each other. 



* Numerous additional multiparameter studies were submit- 
ted after the present paper, the most similar in focus being 
those by the 2dF redshift survey team | p7|j28) . 
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A. CMB data 
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1. CMB data used in our analysis. Error bars do not 
calibration or beam errors which allow substantial vertical 
and tilting for some experiments. 
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FIG. 2. Combination of all data from Figure |l| These error 
bars include the effects of beam and calibration uncertainties, which 
cause long-range correlations of order 5%-10% over the peaks. In 
addition, points tend to be anti-correlated with their nearest neigh- 
bors, typically at the level of a few percent. The horizontal bars 
give the characteristic widths of the window functions (see text). 



Figure [l] shows the 105 band power measurements used 
in our analysis. Starting with the data tabulated in 
we have added the new measurements from CBI [26|, 
QMASK ||, BOOMERaNG @, DASI § and Max- 
ima ||. Since QMASK combines the Saskatoon |^] and 
QMAP |34|-|36fl datasets, these have been omitted. A re- 
cent data review is given in | p7| . 

The success of experimental CMB work has made data 



plots such as Figure [j] increasingly bewildering and dif- 
ficult to interpret. Not only do many data points with 
widely different error bars overlap, but important cor- 
relations due to calibration and beam uncertainties are 
difficult to visualize graphically and tend not to be in- 
cluded in the plotted error bars. The obvious solution to 
this problem is some sort of data compression. 

A radical but common example of this is to simply 
throw away most of the data and show/analyze only one 
or two experiments, often the most recent. This is not 
ideal, however, since it both wasteful of information and 
lacks an objective criterion for data culling. Moreover, 
among the most accurate and thoroughly systematics- 
tested measurements on large scales still come from older 
maps (COBE DMR f| and QMASK ||). 

A more desirable alternative is to average the data 
together somehow, into a single measurement of power 
on each angular scale. Such CMB data c omp ression has 
been performed by many authors, e.g., 10 3l|^0| -^2|, and 
can retain all cosmological information provided that the 
new power bins are narrow enough to not smooth out 
important power spectrum features. However, such com- 
pression throws away any evidence for discrepancies be- 
tween experiments that may have been present in the full 
data set, so it is important to complement the averaging 
by consistency checks. 



B. Combining experiments 

Let us group the power measurements at hand (say the 
105 measurements of ST 2 from Figure |l|) into a vector y. 
We will model this as 



y = Wx + n, 



(1) 



where x is a vector containing the true power spectrum 
coefficients 5Tf up to some sufficiently large multipole 
f max , the window function matrix W encodes the an- 
gular sensitivity of the measurements (the rows of W 
sum to unity) and the noise vector n represents all forms 
of measurement error. We model the errors as random 
with zero mean ((n) = 0) and with a covariance matrix 
N = (nn*) that is the sum of four terms, 



N = N( mcas ) + N^ sca1 ^ + ]Nf( ical ) + N^ bc 



(2) 



corresponding to basic measurement errors, source cal- 
ibration errors, instrument calibration errors and beam 
errors, respectively. In general, all of these errors de- 
pend on the actual power spectrum x, cither through 
sample variance 1 43 1 or because calibration and beam er- 
rors are multiplicative rather than additive. Below we 
will make the approximation that the relative errors are 
small (<gC 1). In this limit, N reduces to a known ma- 
trix independent of x. Explicit expressions for the four 
matrices in equation (Q) are given in Appendix A. 

Given W, N and y, it is straightforward to invert equa- 
tion (pi) to compute an estimate of the underlying power 
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spectrum x. This problem is mathematically identical 
to that involved in CMB mapmaking except that 

the matrices involved are small enough to be trivial to 
manipulate numerically. As our estimator of x we use 

x= [W'N^W] -1 W*N~ 1 y, (3) 

which can be shown to be unbiased ((x) = x), to mini- 
mize the rms noise in each power band and, if the noise 
properties are Gaussian, to retain all information about 
the true power spectrum x from the original data y 
p4j| . The corresponding covariance matrix of the noise 
e = x — x is 

S = (ee*) = [W'N^W]" 1 . (4) 

The resulting power spectrum x is shown in Figure |^ and 
listed in Table 1. The corresponding covariance matrix 
£ is available at 

www.hep.upenn.edu/ ~ max / cmb / experiments .html . 

When computing this spectrum, we did not treat the 
power as an independent parameter at each multipolc. 
Rather, we treated the power spectrum as piecewise con- 
stant, parametrized by its height x\ in the 24 bands listed 
in Table 1. Since our compressed band powers x are sim- 
ply linear combinations of the original measurements y, 
we are able to compute their window functions exactly 
by taking the same linear combinations of the rows of W 
from equation (|]). Our compressed data set can therefore 
be analyzed ignoring the details of how it was constructed 
(ignoring the first column of Table 1), viewed as simply 
the window matrix times the true power. 

The characteristic widths of these window functions 
are reflected by the horizontal bars in Figure || and listed 
in the Table; the exact windows are available on the 
above-mentioned web sitefjj. This correlation matrix in- 
cludes the residual effects of calibration uncertainty and 
beam errors. These long-range correlations are found to 
be moderate, typically of order 5-10% over the peaks, 
which shows that the experiments have to a certain ex- 
tent been calibrated off of each other. In addition, there 
tends to be a slight anti-correlation between neighbor- 
ing points, typically at the level of a few percent, as the 
power spectrum inversion performs a slight deconvolu- 
tion of the input window functions from the experiments 
used. 



T The horizontal location of a data point in Figure |2| corre- 
sponds to the median (50% quantile) of the absolute value of 
its window function. We use absolute values to be pedantic, 
since some windows go slightly negative in places, although 
this makes a negligible difference for the plot. The horizon- 
tal bars to the left and right extend to the 20% and 80% 
quantiles, respectively, indicating the characteristic window 
function width. These quantiles correspond to the full-width 
half-max (FWHM) for a Gaussian window. 



Table 1 - Band powers combining the information from all 105 
CMB data points from Figure [l]. The 2nd column gives the medians 
and characteristic widths of the window functions as detailed in the 
text. The spectrum was computed treating ST 2 as constant in the 
bands listed in the first column. The error bars in the 3rd column 
include the effects of calibration and beam uncertainty. The full 
24 X 24 correlation matrix and 24 X 2000 window matrix are available 
at www .hep.upenn.edu/ ~ max / cmb / experiments .html . 
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One interesting feature of Figure g is that it shows both 
the first and second peak[j] somewhat lower than a large 
fraction of the data. Indeed, the preferred recalibrations 
for all five multiband experiments flagged in the follow- 
ing subsection are downward. To understand the origin 
of this effect, we performed a series of tests where the 
combined spectrum was recomputed with one or more 
experiments omitted. The explanation centers around 
the BOOM98 data, which combines a sharp constraint 
on the relative heights of the first two peaks (even tak- 
ing the beam uncertainty into account) with a relatively 
large overall calibration uncertainty. Since both Maxima 
and DASI (with one exception) have points with small 
error bars below BOOM98 around the second peak, they 
pull the BOOM98 calibration down. QMASK also pulls 
BOOM98 down because of its statistical weight around 
t ~ 100. Although Figure |l| suggests that influence from, 
e.g., TOCO around the first peak might raise BOOM98, 



'One of the Argentinian authors feels that there is still no 
conclusive evidence for multiple peaks in the power spectrum, 
especially in regard to any bets or wagers that may or may 
not be outstanding. 
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this pull is weaker because of TOCO calibration uncer- 
tainties and since error bars are overall smaller out at the 
second peak. This somewhat low normalization persists 
even if any one of Maxima, DASI or QMASK is excluded 
from the analysis. 



C. Comparing experiments 

Since the combined power spectrum presented above is 
only as reliable as the data that went into it, let us now 
test this data for internal consistency. 

As mentioned above, equation (|lj) is analogous to the 
CMB mapmaking problem, which means that all meth- 
ods developed for comparing and combining maps can 
be applied to comparing and combining power spectra as 
well. Given two power spectrum measurements yi and 
y2 modeled as in equation (nh, e.g., as 



y a = W,x - 



y 2 = W 2 x + n 2 , 



(5) 



we wish to know whether they are consistent or display 
evidence of systematic errors. Specifically, is there some 
underlying power spectrum x such that the data sets yi 
and y 2 are consistent with equation (^|)? Let us consider 
the simple case where the two experiments have identical 
window functions, that is, Wi = W 2 . The general case 
can be reduced to this one: In practice, we start by re- 
ducing all experiments to the simple form W = I using 
the deconvolution method described in Appendix D of 

©■ 

Consider two hypotheses: 

Ho: The null hypothesis Ho that there are no systematic 
errors, so that the difference spectrum z = yi — y 2 
consists of pure noise with zero mean and covari- 
ance matrix (zz*) = N = Ni + N 2 . 

Hi : The alternative hypothesis that the difference spec- 
trum z has the same covariance matrix N but a 
non-zero mean m. 



A straightforward variation of the derivation in | 
that the "null-buster" statistic 46 



shows 
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rules out the null hypothesis Ho with the largest average 
significance {v) if H\ is true, and can be interpreted as the 
number of "sigmas" at which Ho is ruled out [Q. The 
case derived in differed in that the mean vanished 
under H\ but that the covariance matrix contained extra 
signal S — the result was of the same form as above, 
but with Q = S. Note that for the special case Q oc N, 
it reduces to simply v = (x 2 — n)/\/2n, where x 2 = 
z'N _1 z is the standard chi-squared statistic. The null- 
buster test can therefore be viewed as a generalized % 2 - 
test which places more weight on those particular modes 



where the expected signal-to-noise is high. It has proven 
successful comparing both microwave background maps 
||,|4] flJII and galaxy distributions [|7|@. Tips for 
rapid implementation in practice are given in |j2| . 

Equation (||) shows that in our case, all weight is placed 
on a single mode m. More generally, the test pays the 
greatest attention to those eigenvectors of Q whose eigen- 
values are large. Consider first the case of calibration 
errors. Then the two measured power spectra are gener- 
ically expected to have the same shape but different nor- 
malizations, so that the vectors (yi) and (y 2 ) are par- 
allel but with different lengths. In other words, S = 
and m = (z) oc Wx, so the mode m that we want our 
test to be sensitive to is shaped like the the expected sky 
power spectrum itself — we make this choice below. Sim- 
ilarly, beam errors show up in a different mode, which (as 
discussed in Appendix A) is shaped like the sky power 
spectrum times I 2 to first order. 
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FIG. 3. Each curve shows the number of standard deviations 
("sigmas") at which a given experiment is inconsistent with all 
others when its power spectrum 8T is multiplied by a constant r. 



D. Test results 

As emphasized by, e.g., Press p9[ , it is important to 
check any data collection for statistically significant out- 
liers. We do this separately for each of the 24 experi- 
ments from Figure [l] as follows. We form the difference 
spectrum 



yi - ryi, 



(7) 



where yi and y 2 are the deconvolved power spectra pro- 
duced from the experiment under study and from all 
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other experiments, respectively, and vary the normaliza- 
tion parameter r. We take calibration and beam errors 
into account in computing the "all other" spectrum, but 
not for the experiment under consideration. The result- 
ing significance level v at which the difference spectrum 
is inconsistent with zero is plotted as a function of r in 
Figure ||. 

All experiments taken together have now detected 
CMB fluctuations at about the 270a level, and the fact 
that all curves except the DASI one asymptote to just 
under that value as r — ► shows that no single experi- 
ment dominates all others in statistical weight. To the 
right, as r — > oo, each curve asymptotes to the signifi- 
cance level at which the experiment in question detects 
signal. If the experiment under study has no calibration 
or beam errors and everything is consistent, v should be 
near zero for r = 1, where it has vanishing mean and 
unit variance ({v) = 0, Av = 1). Only five of the 24 
experiments show a significant difference at the 2cr-level 
(f(l) > 2). Their ^(r)-curves are plotted in Figure [|, 
are seen to be perfectly consistent as well — the relative 
calibration r simply has to be shifted to a different value, 
which is in all cases lower, for which v drops below 2. 

Above we saw that the experiments could be effectively 
recalibrated off of each other. We also computed v for 
the latest BOOMERaNG data on a 2-dimensional grid, 
varying both the calibration and the beam size. How- 
ever, beam information from inter-experiment compari- 
son is substantially less useful at the present time than 
extracting the corresponding calibration information — 
we found the BOOMERaNG beam constraints obtained 
in this way to be weaker than those measured from the 
instrument directly. 

In conclusion, the quantitative tests described above 
show no evidence of inconsistency between the 24 CMB 
experiments when beam and calibration uncertainties are 
taken into account, and the power spectrum shown in 
Figure ^| is consistent with all of them. The dip around 
I = 50, caused mainly by Python V, may warrant further 
investigation to study if a smoother underlying spectrum 
can be consistent with all the data in that range. 



III. IS THE COSMOLOGY STORY 
CONSISTENT? COMPARING AND COMBINING 
DIFFERENT COSMOLOGICAL DATASETS 

In this section we confront the CMB data with other 
cosmological observations, with the goal being both ac- 
curate constraints on cosmological parameters and vari- 
ous cross-checks on the underlying physical assumptions. 
We first map out the subset of the 11-dimensional cosmo- 
logical parameter space from |30| that is consistent with 
CMB, large scale structure (LSS) and Lyman Alpha For- 
est (LyaF) power spectra, with Big Bang nucleosynthesis 
(BBN) and with direct Hubble constant determinations, 
included and excluded in various combinations. We pay 
particular attention to whether these priors are mutually 



consistent or pull in different directions, both in terms of 
which parts of parameter space they pull towards and in 
terms of how much they increase the overall x 2 ■ We will 



return to the consistency issue in Section [V, comparing 
our "concordance model" with various other cosmological 
constraints. 



A. Analysis method 

We employ the multiparameter analysis method de- 
scribed in jS(| and |H| with the following modifications: 
optional inclusion of Lya power spectra, optional dis- 
cretization of the gravity wave contribution to allow ex- 
plicit limits on this parameter, inclusion of CMB beam 
uncertainties as described in Appendix A and refined pa- 
rameter grids to reflect the improved data accuracy. It 
consists of the following steps: 

1. Compute CMB, LSS and LyaF power spectra for 
a grid of models in our 11-dimensional parameter 
space. 

2. Compute a likelihood for each model that quantifies 
how well it fits the data. 

3. Perform 11-dimensional interp- 
olation and marginalization to obtain constraints 
on individual parameters and parameter pairs. 



Our 11 parameters are 



(t, fi k , f2 A , Wd m , Lu h , f v , n s ,n u A s , r, b). 



(8) 



These are the reionization optical depth t, the primordial 
amplitudes A s , rA s and tilts n s , fit of scalar and tensor 
fluctuations, a bias parameter b defined as the ratio be- 
tween rms galaxy fluctuations and rms matter fluctua- 
tions on large scales, and five parameters specifying the 
cosmic matter budget. The various contributions f2j to 
critical density are for curvature fik, vacuum energy 
cold dark matter f2 c dm, hot dark matter (neutrinos) f2„ 
and baryons fib- The quantities u>h = h 2 Qh and u>dm = 
h 2 Vldn\ correspond to the physical densities of baryons 
and total (cold + hot) dark matter (fidm = fi c dm + Qv), 
and f u = fl u /fldm is the fraction of the dark matter that 
is hot. We assume that the galaxy bias b is constant on 
large scales [ |5l| but make no assumptions about its value, 
and therefore marginalize (minimize) over this parame- 
ter before quoting constraints on the other ten. In the 
adaptive mesh spirit, we iteratively refined our parame- 
ter grid to adequately resolve the peak of the likelihood 
function. Our final parameter grids were as follows: 

• r = 0,0.05,0.1,0.2,0.3,0.5,0.8 

• Q A = 0,0.1, ....,1.0 

• f2 k = ±1.0, ±0.8, ±0.6, ±0.5, ±0.4, ±0.3, ±0.2, 
±0.1, ±0.03, 0, truncated so that 

O m = 1 - fi k - fi A € [0.05, 1] 
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• w dm = .02, .05, .08, .11, .13, .16, .20, .50 

• u> h = .003, .015, .018, .020, .022, .025, .03, .04, .07 

• /„ = 0, 0.05, 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1.0 

• n s = 0.5,0.7,0.8,0.9,1.0,1.1,1.2,1.4,1.7 

• nt = _l.o, -0.7, -0.4, -0.2, -0.1, 

• r = 0, 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1.0, 1.4, 1.8, 2.5 

• A s is not discretized 

• b is not discretized 

The parameter r was only discretized when computing 
constraints involving r — it was treated as continuous 
and marginalized over analytically as in |50| when con- 
straining other parameters. We used the original CMB 
data set rather than the compressed one to compute the 
CMB likelihood, since step 2 above is the fastest of the 
three in our analysis pipeline anyway. 

B. Non-CMB data used 

As our LSS data, we follow [|5(J in using the power 
spectrum measured from by the IRAS PSCz survey ]52] ] 
by ]53|| , discarding all measurements on scales smaller 
than 20/i _1 Mpc (fc > 0.3/i/Mpc) to be conservative. 

As our LyoF data, we use the 13 recent power spec- 
trum measurements of p3} ], with an additional 27% cali- 
bration uncertainty common to all points (see also [ j55f ) . 
We compute theoretical predictions for these 13 numbers 
corresponding to each of the ~ 300 million models in our 
database by first computing the matter power spectrum 
P(k) as in [Q, then shifting it vertically and horizon- 
tally (on a log-log scale) to account for the fluctuation 
growth and Hubble parameter shift between z ~ 2.72 
and the present epoch. We used the growA package Q 
for computing the growth factors. 

We quantitatively explore how our results are affected 
by adding various other constraints ("priors"). For BBN, 
we test the prior tub = 0.02 for simplicity, since the error 
bars on the most recent BBN estimates uib = 0.02 ± 
0.002 (95%) from [[[tJ are smaller than our grid spacing. 
For the Hubble parameter, we test the prior h = 0.72 ± 
0.08 from the HST Hubble Key Project §§], assuming 
a Gaussian distribution for simplicity. We also try the 
priors r = 0, f^k = 0, /„ = 0, r = and n s = 1 in various 
combinations. 



C. Basic results 

Our constraints on individual cosmological parameters 
are listed in Table 2 for four cases. Constraints are plot- 
ted in figures | and | for cases 2 and 3. All tabulated and 



plotted bounds are 95% confidence limits^. The first case 
uses constraints from CMB alone, which are still rather 
weak because of a one-dimensional degeneracy coupling 
curvature, baryons, tilt, tensors, dark matter and dark 
energy as described in the following subsection. The sec- 
ond case breaks this degeneracy by combining the CMB 
information with the power spectrum measurements from 
PSCz, and is seen to give rather interesting constraints on 
most parameters. The third case adds the prior assump- 
tions that the Hubble parameter is h — 0.72 ± 0.08 [ p8| , 
which tightens up many constraints, in particular that on 
SIa- The fourth case adds the assumption that the neu- 
trino contribution is cosmologically negligible (/„ ~ 0). 
This is equivalent to assuming that there is no strong 
mass-degeneracy between the relevant neutrino families, 
and that the Super-Kamiokande atmospheric neutrino 
data therefore set the scale of the neutrino density to 
be u v < 10" 3 HQ- 

Table 2 — Best fit values and 95% confidence limits on cosmological 
parameters. Zi on is the rcdshift of reionization and to is the present 
age of the Universe. For the numbers above the horizontal line, the 
central values are the ones maximizing the likelihood (the best fit 
model). For the numbers below the horizontal line, the limits were 
computed from moments as described in the text, so the central 
values are means rather than those for the best fit model. For 
instance, the Hubble parameters for the best fit models are h =0.51, 
0.48, 0.64 and 0.64 for the four columns, respectively, which differs 
from the mean values in the table. 
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° Bayesean 95% confidence limits are in general those that 
enclose 95% of the area. In this paper, we make the approx- 
imation that the boundary of the confidence region is that 
where the likelihood has fallen by a factor e~ Ax 1/2 from its 
maximum, where A\ 2 = 4 for 1-dimensional cases (such as 
the numbers in Table 2) and A^ 2 = 6.18 for 2-dimensional 
cases (such as figures ^ and ^). As shown in Appendix A 
of |pH , this approximation becomes exact only for the case 
when the likelihood has a multivariate Gaussian form. We 
make this approximation to be consistent with the multidi- 
mensional marginalization algorithm employed here (and by 
most other authors), which is equivalent to the integration 
technique only for the Gaussian case. 



6 



For the first 7 parameters listed in Table 2, the 
numbers were computed from the corresponding 1- 
dimensional likelihood functions (plotted in Figure ^ and 
Figure || for the second and third cases). The best fit 
value corresponds to the peak in the likelihood function 
and the 95% limits correspond to where the likelihood 
function drops below the dashed line at e~ 2 of the peak 
value. For the remaining parameters listed, which are 
not fundamental parameters in our 11-dimensional grid, 
the numbers were computed as in []Tof by calculating the 
likelihood-weighted means and standard deviations over 
the multidimensional parameter space. Here the tabu- 
lated limits arc the mean ±2er. 
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FIG. 4. Constraints on individual parameters using only CMB 
and LSS information. The quoted 95% confidence limits are where 
each curve drops below the dashed line. 
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FIG. 5. Like the previous figure, but adding the prior 
h = 0.72 ±0.08. 
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FIG. 6. Constraints in the o;b)-plane. The large 

pink/grey region is ruled out by CMB data alone at 95% confi- 
dence. The two smaller shaded regions become excluded when im- 
posing additional constraints from Hubble constant measurements 
(h = 0.72 ± 0.08) and PSCz galaxy clustering The yellow/light 
grey band shows the BBN constraints from ]57j. The model best 
fitting the CMB+PSCz+/i constraints has x 2 ~ 126 for 127 degrees 
of freedom. 



D. Matter budget 

We will now investigate the parameter constraints in 
more detail, exploring which conclusions come from what 
assumptions. This subsection is centered around the cos- 
mic matter budget (the densities of baryons, cold dark 
matter, hot dark matter, dark energy and curvature) — 
we return to the inflationary parameters n s , n t and r in 
the next subsection. 
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FIG. 7. Constraints in the (n s , cJb)-plane. The yellow/light 
gray band shows the BBN constraints from J>7| . 
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FIG. 8. Constraints in the (u^dmi /i/)-plane. The four curves 
show contours of constant neutrino mass sum. 

1. Constraints from CMB alone 

Figures [| [?] and || summarize some of the key 2- 
dimensional constraints on the matter budget. It is 
seen that the CMB data alone are now for the first 
time (compare, eg, p0| ) powerful enough to close off 
bounded regions in these planes. The low second acous- 
tic peak inferred from the first Antarctic BOOMERaNG 
results |nj was still consistent with a purely baryonic 
Universe, prompting speculation |6^] that an alternative 
theory of gravity might be able to eliminate the need 
to for dark matter altogether. Since the more accurate 
measurements of the second peak height from the new 
BOOMERaNG, DASI and Maxima data give a higher 
value, even the CMB alone now requires a non-zero 
amount of dark matter, at least at modest significance. 

The increased second peak height is also seen to resolve 
a second hot discussion topic of the past year: the CMB 
lower limit on the baryon density has now dropped down 
in beautiful agreement with the BBN prediction. 




2 510 40 100 200 400 600 800 1000 1200 1400 

Multipole 1 

FIG. 9. The main remaining CMB degeneracy. The three solid 
curves show the best fit models for CMB alone subject to the con- 
straint W|, = 0.02 (red), 04, = 0.03 (green) and ajv, = 0.04 (blue), 
spanning the "banana" allowed by the CMB in Figure ^. Dashed 
curves show the scalar and tensor contributions. All models have 
t — fiy — 0. The parameters p = (Q^, Q.^, ^dmi w bi "s, nt) are 
(-0.043, 0.461, 0.116, 0.020, 0.924, 0.982), 
(-0.134, 0.799, 0.050, 0.030, 1.200, 0.884) and 
(-0.160,0.840,0.043,0.040, 1.401,0.955), respectively. 

A third noteworthy result is that the allowed region 
in our 11-dimensional space is shaped like a long and 
skinny rather one-dimensional tube. This is seen clearly 
in the projection onto the (n s , o>b) plane in Figure |?]. The 
physics underlying this CMB degeneracy is illustrated in 
Figure |[ Starting in the generally favored (white) region 
of Figure ^ and moving up to the right, the Universe 
becomes more and more closed (f2k decreases to nega- 
tive values), which would on its own shift the acoustic 
peaks too far to the left. This is compensated by reduc- 
ing the density of dark matter Wdm and increasing the 
dark energy density fl\ and the power spectrum tilt n s 
so that the peak location stays essentially fixed. Through 
a rather spurious coincidence, the resulting changes in the 
various peak heights relative to the COBE scale can be 
almost perfectly reversed by increasing the baryon den- 
sity and adding substantial amounts of gravity waves. 
This is in essence the familiar (fi^, ^a) angular distance 
degeneracy described in many parameter forecasting pa- 
pers p3| , |64| |, with the extra twist that changing toh, i^dm, 
n s and r as well helps minimize the change in the late 
integrated Sachs- Wolfe effect on COBE scales. 

2. Breaking the CMB degeneracy 

Since this degeneracy involves so many parameters, 
prior constraints on any one of them will help break it. 
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In particular, the reason that this degeneracy is not as 
prominent in the recent analyses by the BOOMERaNG 
, DASI p{| and Maxima ||(| teams is because they all 
assumed negligible gravity waves, r = 0. 
Since the Hubble parameter is given by 



h = 



1 — Qjj — f2yV 



(9) 



it decreases sharply as one moves along the degeneracy 
track, dropping as low as 0.3 at the upper right endpoint 
in Figure [?]. Imposing the prior h — 0.72 ± 0.08 therefore 
shrinks the allowed 11-dimensional region substantially, 
among other things tightening the lower limit on dark 
matter in Figure |6| and the upper limit on baryons in 
Figure [?]. 

Because of its sensitivity to n s and u>dm in particular, 
our galaxy clustering data (the PSCz power spectrum) is 
seen to break this CMB degeneracy completely, shrinking 
the allowed CMB "bananas" to almost round regions in 
Figure || and Figure |7[ The effect of adding PSCz is 
particularly dramatic in Figure ||, since CMB alone has 
almost no sensitivity to the neutrino fraction, whereas 
increasing /„ provides a strong suppression of the galaxy 
power spectrum on small scales. 

The model best fitting the CMB+PSCz+/i constraints 
has x 2 ~ 126 for 105 + 32 + 1 - 11 = 127 degrees of free- 
dom. The effective number of degrees of freedom might 
be a few larger than this, since some of the 11 parame- 
ters had little effect, but even if we ignore this, all fits 
are good in the sense of giving reduced x 2_v alues of order 
unity. 

Adding our additional priors uJh = 0.02, r = 0, £lk = 
and r = cause little further change, since they are all 
consistent with the favored results and there are no major 
degeneracies left to break. 

Adding our LyoF data produced effects quite simi- 
lar to those of adding PSCz: breaking the CMB degen- 
eracy and favoring a flat, roughly scale- invariant Uni- 
verse. We have chosen to highlight the effects of the 
PSCz data rather those of the LyoF data in the fig- 
ures since we found its overall constraining power on 
parameters to be slightly stronger. For instance, the 
CMB+LyoF constraints on tilt, curvature and baryon 
density are 0.91±g;^, -0.03±£!g and 0.02lt°;°^, respec- 
tively. Overall, the PSCz and LyoF were found to be 
strikingly consistent in pulling the CMB in the same 
direction, as was also found in fli] ]. Augmenting the 
CMB+PSCz data with the 13 LyoF points increased the 
X 2 for the best fit model by only 10, and the good agree- 
ment between the galaxy and LyaF can also be seen vi- 
sually, in Figure | 




k [l/h- 1 Mpc] 

FIG. 10. The model best fitting the CMB, LSS and LyaF data 
with h and /„ priors is shown together with the two data sets. The 
red and blue curves are the model predictions for the LSS and LyoF 
data, respectively. The blue curve is simply the red one shifted 
vertically and horizontally to account for the fluctuation growth 
and Hubble parameter shift between z ~ 2.72 and the present 
epoch. The overall calibration uncertainty in the LyoF data is 
seen to let the model prediction lie below the data by a constant 
factor. 

We conclude this subsection by summarizing what is 
obtained by adding successively stronger assumptions as 
in Table 2. 

1. CMB alone now gives constraints on most param- 
eters, but they are generally weak because of the 
above-mentioned degeneracy. 

2. Adding the PSCz galaxy clustering constraints 
breaks this degeneracy, resulting in tight con- 
straints on everything except the dark energy den- 
sity Oa- 

3. Adding the constraint h — 72 ±8 finally nails down 
SIa- It also sharpens the limits on fik and f„. 
Adding this one constraint raises x 2 by as much 
as 4 (from 122 to 126), which reflects a slight ten- 
sion between the /i-prior and the CMB+PSCz data, 
which alone prefer the lower range h =0.57^'|J. 

4. Adding an f u — prior boosts x 2 onr y by 0.1, 
further sharpening the f2k, Q\ and r constraints 
slightly, and none of our additional priors (includ- 
ing LyaF) have much of an effect. The only excep- 
tion is imposing n s = 1, which is slightly disfavored 
by the data and raises x 2 by 4. 







E. Inflationary parameters 

We now turn our attention to the parameters associ- 
ated with inflation. Since space remains perfectly consis- 
tent with the inflationary flatness prediction despite the 
sharp error bar reduction caused by the new latest CMB 
data (fi k = 0.00 ±0.06 at 95% for CMB+PSCz+/i), it is 
interesting to quantify the constraints on the parameters 
n g , n t and r to see how they compare with the predic- 
tions of various models, as has previously been done using 
earlier CMB data 0|l4|j27 




0.6 0.8 



Scalar spectral index n s 



FIG. 11. Constraints in the (n s , r) plane. The large green/light 
grey region shows the joint constraints from CMB and LSS. The 
blue/grey region shows the effect of adding the prior h = 0.72±0.08. 
The small yellow/light grey region shows the effect of adding 
f v = 0. The red/dark grey region shows the effect of adding 
the BBN prior = 0.02. The straight line shows the prediction 
r = (200/9)(l — n a ) from a power law class of inflation models. 
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FIG. 12. Same as previous figure but for the (nt,r) plane. 
The straight line shows the "inflationary consistency relation" 
r = -(200/9)n r . 



As discussed in |6q| , there has been a fair amount of 
notational confusion in the literature surrounding the 
tensor-to-scalar ratio r. There are two logical ways to 
define this ratio: either in terms of the fundamental pa- 
rameters of the power spectrum (or, equivalently, of the 
inflationary model space), or in terms of the observables, 
usually the CMB quadrupoles. We adopt the former ap- 
proach, and define 



At 
A' 



(10) 



where A s and A t are the scalar and tensor fluctuation 
amplitudes as defined in |65f| . For inflation models where 
the slow-roll approximation is valid, this ratio is related 
to the tensor tilt n t by the so-called inflationary consis- 
tency condition p5|j66] | 

200 , , 

r= — —n t . (11) 

A power law class of inflation models (see |ll| for a re- 
view) make the additional prediction that 

r=?f(l-n s ). (12) 

i.e., that n s = n t + l. (Although the quantity n t + l would 
be a more natural definition for the tensor spectral index, 
we will stay true to the astronomical tradition of clinging 
to silly notation for historical reasons.) 

A common alternative definition of the tensor-to-scalar 
ratio is the quadrupole ratio 



R 



^jscalar 



(13) 



In this case, the inflationary consistency condition is |6q ] 
R re -6.93n t (14) 



for the special case where f2k = = 0. 
Writing the relation between R and r as 

r ss aR, 



(15) 



the proportionality constant a is typically between 2 and 
5 — it depends on the values of Oa and ilk via the late 
integrated Sachs- Wolfe effect. 

As shown in Table 2, Figure || and Figure ||, the data 
prefers no tensor contribution at all, placing a 95% upper 
limit r < 0.5 for the CMB+PSCz+/i case. For compar- 
ison with prior work, this corresponds to a quadrupole 
ratio of 0.2 in the sense that this is the quadrupole ratio 
for the best fit model in our grid with this r-value (which 
is by definition ruled out at exactly the 95% level). For 
comparison, previous m ultip arameter analyses incorpo- 
rating gravity waves jl3],[l4j found r < 0.3 using older 
CMB data and stronger priors. 
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Figures |ll| and [T^ show the joint constraints on r with 
the scalar and tensor tilts, respectively, and allows com- 
parison with the predictions of equations (11) and (12). 
The constraints in the (n s , r) have becomes progressively 
sharper during the past year Ju],[l4|] , and are now gradu- 
ally becoming quite interesting. In particular, Figure [ll| 
shows that the preference for a slight red-tilt (n s ~ 0.9) 
and low r in the data favors so-called "small-field" models 

& ra 

Figure |12j shows that the holy grail of testing the in- 
flationary consistency relation (^Tj) is still a ways off. As 
expected, the constraints on the tensor tilt n t are stronger 
for high tensor normalizations r and vanish completely 
when r does. 



Consensus has yet to be reached on the second deci- 
mal of the BBN predictions, with the value LOb = 0.023 
from a recent VLT deuterium study [6^ ] lying above 
the 95% range of [p7| . As a reality check, our baryon 
value also agrees with slightly less accurate estimates of 
the abundance in the local Universe — e.g., the range 
0.007 < 51b ^ 0.041 inferred from a low-redshift inven- 
tory Q and the range 0.015 < uib < 0.03 at redshifts 
of a few from the Lya forest [^9|,^0|, although the lat- 
ter tend to be on the high side. Similarly, the inferred 
baryon fraction Wb/wdm ~ 18% agrees with that inferred 
from galaxy clusters p2|]7l| ] , although this match weakly 
prefers lower /i-values. 



IV. DISCUSSION 

In Section |n], we compared and combined the different 
CMB experiments, finding that a consistent picture of 
the angular power spectrum emerged. In Section III, we 



compared and combined a limited number of cosmologi- 
cal data sets (power spectra from CMB, LSS, LyaF and 
various priors), finding that a consistent set of cosmolog- 
ical parameters emerged that provided a good fit to this 
data. 

We conclude with some remarks on how these parame- 
ter measurements match up with the many other cosmo- 
logical observations that probe these parameters, focus- 
ing on the cosmic matter budget. In light of the check- 
ered history of many cosmological parameters, where tiny 
quoted error bars have repeatedly masked larger uncer- 
tainties about underlying assumptions, such end-to-end 
consistency checks are crucial. For instance, the Hub- 
ble parameter h has dropped by a factor of eight since 
Hubble's original paper, the BBN baryon density Co>b has 
risen by 50% in less than a decade, and the favored value 
of the cosmological constant 57a has fluctuated wildly at 
the hands on both theorists and observers. 



A. Baryon density Wb 

The rise of the second peak in the new data has com- 
pletely eliminated the tension between BBN and CMB, 
and they two are now in beautiful agreement that the 
baryon density lj^ ~ 0.02. This agreement was noted 
in the latest team papers as well 
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That one 

method involving nuclear physics when the Universe was 
seconds old and another involving plasma physics more 
than 100,000 years later give the same answers, despite 
involving completely different systematics, is a landmark 
achievement for cosmology. It greatly boosts the credibil- 
ity of the basic cosmological storyline since the Universe 
was a split second old. This sudden agreement is all the 
more impressive given the lack thereof during the past 
year and the ado this generated. 



B. Dark matter and dark energy density 

By now, there are a large number of independent meth- 
ods for probing the dark matter density u>dm or the (al- 
most identical) total density, including studies of the clus- 
ter abundance at various redshifts |fn]-|75||, mass-to-light 
ratio techniques |76| , |7"7j ] , the baryon fraction fib/^m from 
cluster studies |78jj, cosmic velocity fields and redshift 
space distortions. Most of these probes do not measure 
W(j m directly, but fldm — ^dm/h 2 or some combination 
of this with the fluctuation normalization parameter ag- 
Although each comes with potential systematic errors of 
its own, and there is still some internal controversy in the 
velocity and cluster areas Jr9"[ , they are all broadly con- 
sistent with the range cudm = 0.09 — 0.18 listed in Table 
2. Indeed, it is interesting that some of the last evidence 
supporting f2 ro ~ 1 came from velocity fields and redshift 
distortion studies, and that improved data now gives val- 
ues as low as in Table 2 from both velocity fields [{30 81 
and redshift space distortions [ p3p^j83] ] . 

With CMB strongly favoring an essentially flat Uni- 
verse, the low matter density automatically implies a 
high dark energy density. As has been discussed in nu- 
merous papers (e.g., 0Jl^,0), there are now indepen- 
dent ways of reaching this conclusion, making it robust 
to throwing away information on either supernovae la 
g or large-scale structure. The BOOMERaNG team 
argue that even CMB and Hubble constant measure- 
ments now favor Q\ > — we have seen that this conclu- 
sion hinges on the the important additional assumption 
of negligible gravity waves. 



C. Neutrino density 

For the neutrino density w„, we are still a far cry from 
the grand goal of a precision cross-check between cosmo- 
logical and laboratory measurements, since the two have 
so far provided only upper and lower limits, respectively. 
However, the two limits are steadily creeping closer. At- 
mospheric neutrino oscillations J5?| show that there is at 
least one neutrino (presumably a linear combination of 




11 



and v T ) whose mass exceeds a lower limit somewhere 



in the range 0.04-0.08 eV 



1 



Since 



v 

94eV ^ ' 



(16) 



where m, is the mass of the «th neutrino, this corresponds 
to a lower limit w v > 0.0004 - 0.0008, or f v > 0.003 - 
0.01, Our constraints (see Figure ||) give J2 m i < 4.2eV, 
further sharpening the 5.5 eV limit from a careful analysis 
of previous cosmological data Q (see also Q). The 
mass of the heaviest neutrino is thus in the range 0.04 — 
4.2 eV. 

Moreover, the mass-splittigs indicated by both solar 
and atmospheric neutrino data are much smaller than 
4.2 eV, suggesting that all three mass eigenstates would 
need to be almost degenerate for neutrinos to weigh in 
near our upper limit. This means that the upper limit on 
the (almost identical) masses of the three neutrino states 
would be 4.2/3 = 1.4 eV. 

Note that if, as current data suggests, the mixing be- 
tween the flavor eigenstates v e , and v T is not small, 
it is inappropriate to identify the three mass-eigenstates 
v\, vi and z/3 with these flavor eigenstates. For instance, 
the heaviest eigenstate ^3 is likely to be almost a 50 — 50 
mixture of v T and v^. The correct way to phrase our 
upper limit is therefore as a 1.4 eV upper limit on the 
the mass of V3. 

Finally, a caveat about non-standard neutrinos is in 
order. To first order, our cosmological constraint probes 
only the mass density of neutrinos. Our conversion of 
this into a limit on the mass sum assumes that the neu- 
trino number density is known and given by the standard 
model freezeout calculation; 112 cm~ 3 . In more general 
scenarios with sterile or otherwise non-standard neutri- 
nos where the freezeout abundance is different, the robust 
conclusion to take away is simply an upper limit on the 
total light neutrino mass density of 8.4 x 10~ 28 kg/m 3 . 

How do our results compare with those of other recent 
analyses? The analysis most similar to ours is that of the 
BOOMERaNG team ||. A detailed numerical compari- 
son of their results with our Table 2 is very encouraging. 
Despite major differences in both analysis technique (pri- 
ors, parameter space, marginalization method, etc.) and 
data used (that analysis was limited to BOOMERAnG 
and DMR data), both the central values and the error 
bars are very similar for most parameters when taking 
into account that they and we quote la and 2a errors, 
respectively. This indicates that what is being measured 
is really borne out loud and clear by the data in a way 
that is robust towards data selection or analysis details. 
Perhaps we are inevitably approaching the dreaded day 
when not only cosmology is consistent, but cosmologists 
are as well. 
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APPENDIX A: LINEARIZED MODELING OF 
CALIBRATION AND BEAM ERRORS 

This appendix describes our modeling of the four terms 
of the error matrix from equation (|^), assuming that all 
relative errors are small (<C 1). 

N-(mcas) reflects the part of the errors which are uncor- 
rected between the different experiments and is due to 
detector noise and sample variance. We approximate it 
by 



N 



(meas) 



(Al) 



where ai is defined as the average of the upper and lower 
error bars quoted for di = ST 2 . 

The last three terms in equation (^) reflect correlations 
between measurements due to calibration and beam er- 



rors. As in 1 31 43], N( lcal ) is the part specific to a single 
multi-band experiment and N^ sca1 ^ is the part that is cor- 
related with other experiments that are calibrated off of 
the same (slightly uncertain) source. TOCO, MSAM, 
CBI and BOOMERaNG all calibrate off of Jupiter. To 
be conservative, we assume that the full 5% calibra- 
tion uncertainty from Jupiter's antenna temperature is 
shared by these experiment. The true correlation should 
be lower, since the four experiments observed Jupiter 
at different frequencies. The remaining multi-band ex- 
periments in Figure [l] should not have any such inter- 
experiment correlations. This contribution to the noise 
matrix is therefore 



(A2) 



where 



f 5% if i and j refer to a Jupiter-calibrated point, 



otherwise. 



(A3) 



The factor of 2 in equation (A2) stems from the fact that 
the percentage error on 5Tf is roughly twice that for 5Ti 
as long as it is small. We simply use the observed values 
for Di in this expression. 

Similarly, the remaining calibration term is 



Ng 5 al) = (2r t 



(A4) 



where = if i and j refer to different experiments. 
If band powers i and j are from the same experiment, 
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then Tij is the quoted calibration error with the source [21] 
contribution Sij subtracted off in quadrature. We use 



10% for QMASK, 14% for Python V, 8% for Viper, 8.7% 
for Toco 97, 6.2% for Toco 98, 10% for QMASK, 6.4% 
for BOOM97, 10% for BOOM98, 4% for Maxima and 4% 
for DASI. 

A Gaussian beam of standard deviation 6 suppresses 
power by a factor well approximated by e^ e ^ . Taylor 
expanding this expression shows that a small beam error 
A9 causes the power to be mis-estimated by a percentage 
2£ 2 6A6. Defining 



.Aft 



(A5) 



(A6) 



we can therefore write the beam error as 

N J cam) = hb 3 

when i and j refer to the same experiment, zero other- 
wise. We use this approximation f or BO OM98, where 
9 = 12.97c, 56 = 1.47c, and c = \/81n2 is the familiar 
FWHM conversion factor. The other experiment report- 
ing important beam uncertainties is Maxima, for which 
we use the ^-coefficients published in || in place of the 
approximation of equation ( A5). 
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